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Bloch oscillations (BOs) of polychromatic beams in circularly-curved optical waveguide arrays are smeared 
out owing to the dependence of the BO spatial period on wavelength. Here it is shown that restoring of the 
self-imaging property of the array and approximate BOs over relatively broad spectral ranges can be achieved 
by insertion of suitable lumped phase slips uniformly applied across the array. 
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Photonic lattices have provided in recent years a lab- 
oratory tool to visualize optical analogues of coherent 
phenomena generally encountered in solid-state physics, 
such as Bloch oscillations (BOs) [T}S] and dynamic lo- 
calization (DL) [5l[6| in lattices driven by external dc 
or ac fields. In waveguide arrays, a common engineer- 
ing approach to realize optical BOs and DL is to tailor 
the local curvature of the array axis [7|l8j. DL and BOs 
have been observed in bent waveguide arrays with either 
periodically- varying or constant curvature to mimic ac or 
dc fields [SliniinillO]. In optics, BOs and DL enable light 
diffraction management, with the possibility of periodic 
image reconstruction [8l[TT]. Such a self-imaging effect, 
however, is generally limited to relatively narrow spec- 
tral beams. In case of DL, destruction of self-imaging 
for broadband beams stems from the resonant nature 
of the DL condition [12 . Here, periodic self-imaging is 
exactly fulfilled solely for a target wavelength, whereas 
other spectral components experience residual diffrac- 
tion. In |13J, Garanovich and coworkers suggested spe- 
cial profiles of waveguide axis curvature that mix first 
and second resonances of DL, leading to approximate 
DL over an extremely broad spectral region. The exper- 
imental demonstration of broadband DL has been very 
recently reported by Szameit and collaborators [H]. In 
case of BOs, destruction of exact self-imaging with poly- 
chromatic beams has a different reason. Here periodic 
beam reconstruction is exactly attained for each spec- 
tral beam component, however the self-imaging period, 
given by the Bloch period sb^ turns out to be propor- 
tional to the wavelength A: therefore, for a beam spec- 
trally broadened by less than ~ 10% around its carrier, 
field reconstruction is fully lost just after few BO cycles. 
In this Letter it is shown that polychromatic optical BOs 
can be approximately achieved over a broad spectral 
range by suitable insertion of lumped phase slips in the 
array, which combat the dispersion of BO period with 
wavelength. Let us consider light propagation in low- 
contrast-index one-dimensional waveguide arrays with a 
weakly curved axis in the (x, z) plane. In the curvilinear 
coordinate system (5,7^) of Fig. 1(a), where s is the arc 
length of the curved waveguide axis, the propagation for 



the spectral field component ijjiri^s^X) at wavelength A 
is described by the Schrodinger-like equation 
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where h = A/(27r) is the reduced wavelength, Ug is the 




Fig. 1. (Color online) (a) Schematic of a curved waveg- 
uide array, (b) Refractive index profile n{r]) — rig of 
the array used in numerical simulations (a = 8 /im, 
Ug = 1.45). (c) Behavior of r(A), entering in Eq.(3), 
versus wavelength for (i) the 8-cm-long circularly-curved 
array without lumped phase slips (solid curve); (ii) ai = 
<^3 = Ao/(2ans) and a2 = (dashed curve); (iii) for 
ai = a2 = as = Xo/{3ans) (dotted curve), (d) Behavior 
of waveguide coupling rate A versus wavelength. 

bulk refractive index, n{r]) is the periodic index profile 
of the array with lattice period a, and i<i{s) is the lo- 
cal curvature of waveguide axis. In the single-band and 
nearest-neighbor tight-binding approximations, Eq.(l) 
reduces to standard coupled mode equations for ampli- 
tudes Cn{s;X) of modes trapped in the n-th waveguide 
of the array [8l[l3] 

A / A \ / \ 27Tnsf^(s)an , . 

ICn = -A(A)(Cn+l + Cn-l) H Cn (2) 

where A (A) is the coupling rate, at wavelength A, be- 
tween adjacent waveguides, and the dot denotes the 
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derivative with respect to the arc length. Here we con- 
sider a waveguide array with constant curvature = 
l/R^ into which a sequence of lumped phase gradients 
are superimposed at the arc lengths 5i, S2,.---, i-e. we as- 
sume in Eq.(2) f<i{s) = 1/ R-\-aiS{s — si)-\-a2S{s — S2)-\-...^ 
where the parameters ai, a2,.-- define the phase slips cpi^ 
(^2v introduced uniformly across the array at the arc 
lengths si, S2,.--- The relation between the phase slip (pi 
and the parameter ai is simply obtained after integration 
of Eq.(2) in the infinitesimal interval {s = = s^), 
yielding cpi = 27rnsaai/ X. In practice, a lumped phase 
gradient can be introduced (i) by a sudden tilt of waveg- 
uide axis [11] by a small angle 6; in this case ai = 6 and 
(pi — tiO/Ob^ where Ob = X/{2ans) is the Bragg angle; 
(ii) by local waveguide segmentation, channel narrowing 
or index change modulation [15] . To investigate the self- 
imaging property of the arrayed structure, it is enough 
to consider the impulse response of the array correspond- 
ing to excitation of the waveguide n = 0. The solution to 
Eq.(2) with the initial condition Cn(0) = (5^,0 reads [12] 

|c„(s)p = J^(2r(A)), (3) 

where Jn is the Bessel function of order n, r(A) = 
A{X)\q{s,X)\ and 

q{s,X) = / (i^exp[-i7(0/A], 7(5) = 27rnsa / d^K.{^). 
Jo Jo 

(4) 

The condition for self-imaging after a propagation length 
5, i.e. |cn(s)p = ^n,o, is thus q{s^X) = 0. For the BO 
problem in absence of phase slips, k{s) = 1/R and 
\q{s^ A) I = s \ sm{7TS / s b) \ / {tts / s b) ^ where 



sb{X) = XR/ijisa) 
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is the BO period for the spectral field component of 
wavelength A. For a monochromatic beam with carrier 
wavelength Aq, periodic self-imaging is thus attained at 
propagation lengths s which are integer multiples of 
the BO cycle sb{Xo). Unfortunately, for spectrally-broad 
beams different wavelength components A show a shifted 
BO cycle 5b(A), making self-imaging and polychromatic 
BOs to rapidly washing out for a relatively broad spec- 
trum of the incoming beam. As an example. Fig. 2 (a) 
shows the breakdown of self-imaging and smearing out of 
a breathing BO mode for polychromatic light in a L = 8- 
cm-long circularly-curved array as obtained by numeri- 
cal simulations of Eq.(l) for an index profile of the array 
shown in Fig. 1(b) and for a radius R = 4Q cm. The inten- 
sity distribution of Fig. 2 (a) is obtained by the incoher- 
ent superposition of five monochromatic fields of equal 
power (at A =540, 560, 580, 600, and 620 nm), which 
excite simultaneously the waveguide n = 0. For com- 
parison. Fig. 2(b) shows the evolution of beam intensity 
for the carrier spectral component solely (A = Aq = 580 
nm), which undergoes periodic self imaging with spatial 
period sb{Xq) = XoR/{nsa) = 2 cm. Note that the full 
propagation length L cm comprises exactly = 4 BO 



cycles at the carrier wavelength. The washing out of self 
imaging at the output plane for a spectrally broad beam 
is related to the rapid deviation from zero of r(A) as A 
varies by a a few percents around the reference wave- 
length Ao, as shown by the solid curve in Fig. 1(c). The 
allowed spectral band AA for (approximate) self-imaging 
can be roughly estimated from the condition 2|r(A)| < 1 
[see Eq.(3)], and turns out to be A A ~ 20 nm for the 
solid curve of Fig. 1(c). The introduction of appropriate 
phase slips across the array may be exploited to flatten 
the function r(A) at around A = Aq. To this aim, let us 
consider {N — 1) lumped phase gradients at si = 5b (Aq), 
S2 = 2sb{Xo),'",sn-i = {N — 1)sb{Xo), i.e. let us assume 
K,{s) = 1/^ + Xl^Jo^ (^k^{s — Sk) with So = and ao = 
for definiteness. From Eq.(4) with s = the expression 
of q can be calculated in closed form and reads 

L , . sin(7rAo/A) 
q{X) = - exp (-*7rAo/A) ^^J^ ' F{7rX/Xo) (6) 

where we have set F{^) = ^^Sq ex.p{—2i^pj.) and pk = 
k-\- {aus/ Xo){ao +ai + ...a/c). Appropriate choice of the 
values ak may flatten the shape of r(A) at around its 
zero A = Ao- A first approach consists of considering 
the Taylor expansion of r(A) at around Ao and imposing 
F = (dV/dX) = ... = {d^T/dX^) = at A = Ao up to a 
certain order M > 1, which defines the spectral flatness 
order. This requires F(7r) = F'{7r) = ... = F^^-^^ir) = 
0, which yields 
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Fig. 2. (Color online) Intensity beam evolution (a) in the 
polychromatic, and (b) monochromatic (at carrier wave- 
length) regimes in a L = 8-cm-long circularly-curved ar- 
ray without lumped phase slips. The insets at the top 
are the intensity profiles at the output plane. Parameter 
values are given in the text. 

for / = 0, 1, .., M — 1. As an example, let us consider 
the case = 4, previously discussed in Fig. 2. Then one 
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Fig. 3. (Color online) Polychromatic BOs as obtained 
by insertion of two lumped waveguide tilts at 5i = 2 cm 
and S3 = 6 cm (tilt angle ^ = = 25 mrad). In (b) 
the output intensity profiles at the output plane for the 
different spectral components are depicted. 
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Fig. 4. (Color online) Same as Fig. 3, but with the in- 
sertion of three lumped waveguide tilts at 5i = 2 cm, 

52 = 4 cm and S3 = 6 cm (tilt angle = 20b — 16.7 
mrad) . 

can check that Eq.(7) is satisfied up to the order M = 2 
by simply assuming ai = as = Xo/{2ans) and a2 = 0. 
This corresponds to the introduction, at the arc lengths 
si = sb(Ao) and S3 = 3sb(Ao), a lumped tt phase shift 
uniformly across the array. The corresponding behavior 
of r(A) is shown by the dashed curve in Fig. 1(c). Poly- 
chromatic BOs and (approximate) self-imaging over the 
entire spectral range (540,620) nm is demonstrated in 
Fig. 3(a), where the intensity evolution of the polychro- 
matic beam of Fig. 2 (a) is now computed with insertion 
of the phase gradients at si = sb(Ao) = 2 cm and at 

53 = 3sb(Ao) = 6 cm. In the numerical simulations, the 
lumped TT phase slips are simply obtained by tilting the 
waveguides (tilting angle = Ob = 25 mrad) as in [TT]. 
The spectrally-resolved output intensity profiles of the 
polychromatic beam are also shown in Fig. 3(b). Note 
that, with this choice of ai, a2 and a3, the function r(A) 
remains reasonably small in the range A ~ (540, 620) 
nm, however it rapidly increases outside this interval. A 



different choice for ai, a2 and a^^ which broadens the 
spectral range where 2r(A) remains smaller than ~ 1, 
is ai = a2 = Oi3 = Ao/(3ans). In this case, r'(Ao) and 
F (Ao) do not vanish, however two additional zeros of 
F are introduced at A :^ 516 nm and A ^i:^ 618 nm, as 
shown by the dotted curve in Fig. 1(c). The correspond- 
ing behavior of polychromatic BOs, for a beam covering a 
spectral extent AA ~ 130 nm, is shown in Fig. 4. Different 
algorithms, based on optimization methods (e.g., trying 
to minimize J (iA|r(A)|), might be investigated to flatten 
the function r(A) over a prescribed spectral range. As a 
general rule, the increase of the number N of phase gra- 
dients makes it the flattening procedure more flexible, 
and thus the polychromatic imaging of higher quality or 
applicable to a broader spectral range (for instance by 
increasing the spectral flattening order M). 
In conclusion, broadband self-imaging based on poly- 
chromatic BOs in circularly-curved waveguide arrays has 
been proposed by insertion of lumped phase slips, which 
combat the dispersion of BO period for the various spec- 
tral beam components. 
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